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^ . Abstract 

We consider a homogenization of elliptic spectral problem stated in a perfo- 
rated domain, Fourier boundary conditions being imposed on the boundary of 
Oh I perforation. The presence of a locally periodic coefficient in the boundary op- 

erator gives rise to the effect of a localization of the eigenfunctions. Moreover, 



< 

' ■^ \ the limit behaviour of the lower part of the spectrum can be described in terms 



(N 



of an auxiliary harmonic oscillator operator. We describe the asymptotics of the 
eigenpairs and derive the estimates for the rate of convergence. 
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try 1 Introduction 



The paper deals with a spectral problem for a second order divergence form 
elliptic operator in a periodically perforated bounded domain in M. d . Assuming 
that on the perforation border a homogeneous Fourier boundary condition is 
stated, and that the coefficient of the boundary operator is a function of "slow" 
argument, we arrive at the following eigenvalue problem 



-div(a(x/e)Vu £ (x)) = X £ u £ (x), x G fi e , 
a(x/e)Vu £ (x) • n = — q(x)u e (x), x 6 £ e , (1-1) 
u £ {x) = 0, xedfl; 

here e is a small positive parameter defined as a microstructure period. 

We impose some natural regularity and connectedness conditions on the per- 
forated domain Q e , as well as usual periodicity and uniform ellipticity conditions 
on the matrix a(y). These conditions are specified in detail in the next section. 

Our crucial assumptions are 
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• q € C 2 (ft), and q(x) > q > in Q. 

• The function q has only one global minimum point in fi. The global mini- 
mum is attained at an interior point of 0. 

• The Hessian matrix d 2 q/dx 2 evaluated at the minimum point is positive 
definite. 

Under the first two assumptions the localization phenomenon holds. Namely, 
for any k £ N the k-th eigenfunction of problem (jl.ip is asymptotically localized, 
as £ —7- 0, in a small neighbourhood of the minimum point. In particular, the 
properly normalized principal eigenfunction converges to a (5-function supported 
at the minimum point. 

In the paper, assuming that all the above conditions are fulfilled, we construct 
the first two leading terms of the asymptotic expansions for the k-th eigenpair, 
k = 1,2,... . 

These asymptotic expansions have a number of interesting features. First of 
all, the mentioned expansions are in integer powers of e 1//4 . Then, the localization 
takes place in the scale e 1//4 . In this scale the leading term of the asymptotic 
expansion for the k-th eigenfunction proved to be the k-th eigenfunction of an 
auxiliary harmonic oscillator operator. 

If q G C 3 (J7), then we also obtain the estimates for the rate of convergence. 

We suppose that q does not oscillate just for presentation simplicity. The 
techniques developed in the paper also apply to the case of locally periodic coef- 
ficients q = q(x,x/e), a = a(x,x/e) with q(x,y) and a(x,y) being periodic in y, 
see Remark 12.11 and Theorem 12.21 below. 

Previously, the localization phenomenon in spectral problems has been ob- 
served in several mathematical works. In [1] the operator with a large locally 
periodic potential has been considered. The localization appeared due to the 
presence of a large factor in the potential and the fact that the operator coeffi- 
cients depend on slow variable. 

In [2] the Dirichlet spectral problem for the Laplacian in a thin 2D strip 
of slowly varying thickness has been studied. Here the localization has been 
observed in the vicinity of the point of maximum thickness. The large parameter 
is the first eigenvalue of ID Laplacian in the cross-section. This eigenvalue grows 
to infinity because the thickness of the strip asymptotically vanishes. 

In the mentioned works, under natural non-degeneracy conditions, the asymp- 
totics of the eigenpairs was described in terms of the spectrum of an appropriate 
harmonic oscillator operator. However, the localization scale was of order y/e 
with e being the microscopic length scale. 

The localization in the scale e 1 / 4 that is observed in the present paper, is 
not standard. It should also be noted that although the operators in (jl.ip do 
not contain a large parameter, such a parameter is presented implicitly because 
(d — l)-dimensional volume of the perforation surface tends to infinity. 
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The homogenization of spectral problem with a constant or periodic 

functions a and q has been addressed in [7j. 

Spectral problems in perforated domains with Dirichlet and Neumann bound- 
ary condition at the perforation border are now well studied. There is a vast 
literature on the topic, see, for instance, [8], [6]. 

In the paper we combine asymptotic expansion techniques with various vari- 
ational and compactness arguments and scaled trace and Poincare type inequal- 
ities. 

2 Problem statement 

We start by describing the geometry of the domain. Let K = [0, l) d and E cM. d 
be a K-periodic, open, connected set with a Lipschitz boundary S; the comple- 
ment W 1 \ E is denoted by B. We also assume that K n E is a connected set, 
and K n B <s K, so that B = M. d \ E consists of disjoint components. In what 
follows, Y = K n E denotes the periodicity cell, and S° = K n dB = K n S the 
boundary of the inclusion. The symbols \Y\d and |£°|<2-i stand for the measures 
of Y and the (d — l)-dimensional surface measure of E°, respectively. 

For every i G Z d we denote Y* = e(i + Y), = eS f| 3^, and B l £ = eBf] Y £ \ 
Given fi, a bounded domain in H.™ with a Lipschitz boundary we introduce 
the perforated domain 

n £ = n\\jBi, i £ = {iez d : y* c n}. 

Notice that 17 £ remains connected, the perforation does not intersect the bound- 
ary <9$7, and 

dn e = dn{Jx £ , s £ =(Js*. 

iSle 

In the perforated domain 17 £ we consider the following spectral problem: 

-div(a £ (x)Vu £ (x)) = X £ u £ (x), x £ Q £ , 

< a £ (x)Vu £ (x) • n = — q(x)u £ (x), x £ S £ , (2-1) 

_ u £ (x) = 0, x G 90. 

Here e is a small positive parameter, a £ {x) = a(x/e) with a(y) being & d x d 
matrix, n is an outward unit normal; the usual scalar product in R d is denoted 
by " •". 

In the sequel we assume that the following conditions hold true: 

(HI) a{y) is a real symmetric d x d matrix satisfying the uniform ellipticity 
condition 

d 

for some A > 0. 
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Figure 1: Domain Q E 

(H2) The coefficients ciij(y) are Y-periodic and, moreover, a,ij(y) G L c 
(H3) The function q(x) G C 3 (R d ) is positive. 

(H4) The function q(x) has a unique global minimum attained at x = G Q. 
Moreover, in the vicinity of x = 

q(x) =q(0) + - x T H{q) x + o(\x\ 2 ), 

with the positive definite Hessian matrix H{q). 
It is convenient to introduce the notation 

H^(Fl e ,dn) = {u£ JT^fie) : u = on dfl}. 

The weak formulation of spectral problem (|2.ip reads: find A e G C (eigenvalues) 
and it e G .Hg(f2 e , 90), u 6 / 0, such that 

J a £ Vu e -Vvdx + J qu e vda = X e J u £ vdx, v G Hq(Q). (2.2) 



Lemma 2.1. For any e > 0, i/ie spectrum of problem (|2.2j) is real and consists 
of a countable set of points 

< X\ < X 6 2 < ■ ■ ■ < X £ < ► +00. 

Every eigenvalue has a finite multiplicity. The corresponding eigenf unctions nor- 
malized by 

u\ u £ j dx = 5{j , 



form an orthonormal basis in L 2 (S\). 

We omit the proof of Lemma 12.11 which is classical. 

Under the assumptions (H1)-(H4) we study the asymptotic behaviour of 
eigenpairs (A e , u £ ), as e — > 0. 

To avoid excessive technicalities for the moment, we state our main result 
in a slightly reduced form, without specifying the rate of convergence. For the 
detailed formulation of the main result see Theorem [ 
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Theorem 2.1. Let conditions (H1)-(H4) be fulfilled. If {X £ -,u £ -) stands for the 
jth eigenpair of problem (12 -If) , then for any j, the following representation takes 
place: 



X s _ 1 |E l 



d-1 



q(0) + ^= £ , u £ (x)=v £ (- TJI ) 



where (/z|, uj(z)) are such that 

• /ij converges, as e — > 0, to the jth eiegnvalue fij of the effective spectral 
problem 

- div(a cff Vt>) + (z T Qz) v = pv, v£ L 2 (R d ), (2.3) 
where a cH is a positive definite matrix (see (|3.19|) ); Q is defined by 

Id— 1 



Q 



1 |E°i 



2 \Y\ d 

with H{q) being the Hessian matrix of q at x = 0. . 

• If Hj is a simple eigenvalue, then, for small e, fJu is also simple, and the 
convergence of the corresponding eigenf unctions ( extended to the whole M. d ) 
holds 

\\v e j -Vj\\ L 2 (Rd) ->0, s->0. 

Remark 2.1. Theorem \2.1\ can be generalized to the case of locally periodic 
coefficients in h2.1\) . 

Namely, let us consider the following problem: 
( 



with 



Assume that 



-div(a £ (x)Vu £ (x)) = X £ u £ (x), x E J2 e , 
a £ (x)Vu £ (x) ■ n = —q £ {x) u £ (x), x 6 T, £ , 
u £ (x) = 0, x G dtt 

a £ {x) = a(x, x/e), q £ (x) = q(x,x/e). 



(2.4) 



Qij(x,y) andq(x,y) are Y -periodic in y functions such that Oy(x,y), q(x,y) £ 



C 2 ' a (R d ;C a (Y)) with some a > 0. 

The matrix a(x, y) satisfies the uniform elipticity condition. 
The local average of q defined by 



q(x) 



1 



|E°| 



d-i 



(x,y) da y , 



admits its global minimum at x = 0. 
In the vicinity of x = 

q(x) = g(0) + \x T H(q) x + o(|a;| 2 ) 
with the positive definite Hessian matrix H(q). 
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• x = is the only global minimum point of q in 0,. 
Then the following convergence result holds. 

Theorem 2.2. //(A|,iip stands for the jth eigenpair of problem (|2.4p . then for 
any j, the following representation takes place: 

where (/j,j,Vj(z)) are such that 

• /ij converges, as e — > 0, to the jth eiegnvalue [Aj of the effective spectral 
problem 

- div(a cff Vt;) + (z T Pz) v = /it), v£ L 2 (R d ), (2.5) 

where 

and a cff is a positive definite matrix defined by 

1 



a cff 



13 \y\ d 

Y 



<*k(0,C)(*jy + 9k#i(OK> 



with the functions Nj solving auxiliary cell problems 

-div c (a(0,C)V c iV fc (C)) = div c a. fc (0,C), k = l,...,d, (£Y, 
a(0,()V ( ;N k -n = -a ik (0,C)n h ( G E°, 

N k (c) g ffi (y), 

7f a simple eigenvalue, then, for small e, fij is also simple, and the 

convergence of the corresponding eigenf unctions ( extended to the whole M. d ) 
holds 

\\ v j ~ "jIIl^r") -> °> e^O. 



3 Proof of Theorem 2.1 



3.1 Preliminaries. Estimates for Af 

In this section we estimate the first eigenvalue Af of problem f|2. 1 1) . To this end we 
use the variational representation for Af . Let us recall that, due to the classical 
min-max principle (see, for example, [3]), 

/ a £ Vv -Vvdx + fq(v) 2 da 
Af= inf Of 5f . (3.1) 



Lemma 3.1. The first eigenvalue of the spectral problem (|2.ip satisfies the esti- 
mate 

- 9(0) + 0(1) < Af < - q(0) + 0( £ -V2) ; e _> 0. 
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Proof. We start by proving the estimate from below. By (|3.ip . 

Af > inf \ [ a £ Vv-Vvdx + q(0) f(v) 2 da\. 
veH^ne.dn) ^ J J > 

The last infimum is attained on the first eigenfunction of the following spectral 
problem 

— div(a £ (x)Vw £ (x)) = v £ w £ (x), x G O e , 
a £ {x)Vw £ {x) ■ n = — q(0)w £ (x), x G S e , 

w £ (x) = 0, x G d£l. 

It has been proven in [7] that the first eigenvalue of this problem admits the 
following asymptotics: 

^i = -^Sp<K0) + o(i), e^o. 



Thus, 



Af >il3p g (o) + o(i), e^O. 

£ •» d 



We proceed to the derivation of the upper bound for Af . Choosing v G Cg°(fl) 
as a test function in (|3.ip . one can obtain a rough estimate 

Af < Ce-\ (3.2) 

with a constant C independent of e. To specify C one should choose a "smarter" 
test function. Let us take v G Co°(R d ), \\v\\u2(md\ = 1, and choose v(x/e a ) as 
a test function in (|3.ip . < a < 1/2. Note that if suppw C Br(0), for some 
R > 0, then suppv(x/e a ) C B £ a R (0). Then we obtain 

/g(a;)|i>(^)| 2 da + 0(e- 2a e da ) 
Af < - ^ • 

/ KlOl dx 

Taking into account assumption (H4) and using Lemma 14.11 one has 
IHYT 1 ?(°) / KfO I 2 dx + O^ 2 "" 1 £dQ ) + 0(e- 2a e da ) 

I Id „ 

Af < ^ ■ - 2 . 

/ \v{§;)\ dx 

Notice that the best estimate is obtained for a = 1/4. Finally, 

Af<-^Pg(0) + 0(^/2), £ ^o. (3.3) 

£ 1 1 \d 

□ 

Remark 3.1. When deriving the upper bound for Af , we used a test function 
which is concentrated at x = 0. Namely, the test function of the form v(e~ 1 ^ i x). 
This observation turns out very helpful for the construction of the asymptotics 
of eigenpairs (X £ ,u £ ). 
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The next definition explains the notion of concentration. 

Definition 3.1. We say that a family {w £ (x)} £> o with < c\ < ^\w^\L 2 (Q. e ) — c 2 
is concentrated at xq, as e — > 0, if for any 7 > there is eq > such that 

J \w e \ 2 dx < 7, for all e G (0,£o)- 
Here B~(xq) is a ball of radius 7 centered at xq. 

Lemma 3.2. The first eigenfunction u\ of problem (|2.1j) is concentrated in the 
sense of Definition \3.1\ at the minimum point of q{x), that is at x = 0. 

Proof. Assume that itf, normalized by H^J-f ||z^ 2 (ri e ) = 1> is n °t concentrated at 
x = 0. Then, there exists 7 > such that, for any £q, we have 

\u\\ 2 dx>~{ (3.4) 
n £ \B 1 (o) 

for some e < Eq. 

Estimate (13. 2p together with (|3.1f) imply the estimate 

J \Vu\\ 2 dx < Ce- 1 . 

Then, using Lemma 14.11 we obtain 
Af = J a £ Vu £ 1 -Vu £ 1 dx + ^ l^" 1 j q\u\\ 2 dx + 0{£- l/2 ) 

1 |EVi . f a 

> — — — mm q / ti, ax 

e \Y\ d n £ \B 7 (o) J 

Q e \ B 7(0) 



1 |S°| d 



-{g(0) I \u\\ 2 dx+ f (q(x) - g(0))|<| 2 dx] + 0{e- 1 ' 2 ) 



e \y\d 

n £ ns 7 (o) Q e ns 7 (o) 
Since x = is the global minimum point of q(x), then 

A e > i^Vir min /" u|i2 da . + (0) /" |^|2 dx | + ( £ -l/2). 

e \Y\ d ln e \B 7 (o) J J J 

n £ VB 7 (o) O e nB 7 (0) 

By (HI, 

Af > ~ 9(0) + J l^jfi ( min q - g(0)) 7 + O^ 1 / 2 ), (3.5) 

e |r| d e \Y\ d Vn E \-B 7 (o) / 

that contradicts (|3.3p . Lemma is proved. □ 
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Remark 3.2. The min-max principle allows us to compare the eigenvalues of 
Dirichlet, Neumann and Fourier spectral problems. Namely, denote by X £ D k the 
kth eigenvalue of the Dirichlet problem (u e = on T, e ), and by X £ N k the kth 
eiegnvalue of the Neumann problem (the case q = in (|2.ip ). Then, one can see 
that 

X £ Nik <X £ k <X £ D ,k, fe = l,2,-.-. (3.6) 

It is well-known (see that X £ N k = 0(1) and X £ D k = 0(s~ 2 ), e — > 0. Lemma [3A\ 
specifies estimate (|3.6p for the first eigenvalue Af . 

3.2 Change of unknowns. Rescaled problem 

For brevity, we denote 

Kid ^ \l \d 

where H (q) is the Hessian matrix of q at x = 0. 

Note that Lemma 1331 suggests to study the asymptotics of (A| — e~ l >c(0)), 
rather than of A| itself. On the other hand, when deriving the upper bound in 
Lemma [3 -H we used the test function ^(x/e 1 / 4 ), which allowed us to get the "op- 
timal" estimate. Bearing in mind these two ideas, we first subtract e _1 x(0) u £ (x) 
from both sides of the equation in (|2.ip . and then make the change of variables 
z = e~ x l^x in (|2.ip . Then, the rescaled problem is stated in the domain 

n £ = e- l l^ £ , £; = e- 1 / 4 S £ , 

and takes the form 

-dw(a £ (z)Vv £ (z)) -^v £ = v £ v £ (x), z G i\, 
a £ {z)Vv £ {z)-n = -e l l A q{e l / A z)v £ {z), z G S^, (3-7) 

k v £ {z) = 0, z£ e" 1 /^. 



Here 



v 



£ (z)=u £ (e l ^z), a £ (z) = a(^ /I ), fI £ = ^(X £ -^). (3.8) 



The weak formulation of problem (|3.ip reads: find (/i £ , v £ ) G M.xHq(Q £ ,e 1//4 (9f2), 
w e / 0, such that 

w £ (v £ ,w) = fi £ {v £ , w ) L2( ^ r Vw g flj(n^,e -1 / 4 an). (3.9) 

Here the bilinear form W £ (u, v) is given by 

W £ (u,v) = J a £ Vu-Vvdz - J uvdz + e 1/4 j q(e l/i z)uv da z . (3.10) 
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Remark 3.3. About the extension operator For all sufficiently small e, there 
exists an extension operator 

P £ : Hl&ce-^dQ) -> H^{e- 1/4 n) 

such that 

\\ pe v\\ L 2 {e -i/± n) <C\H L 2 { n e y \\V{P £ v)\\ L 2 {£ -i/4 Q) <C\\Vv\\ L2{ ^ r 

where C is a constant independent of e. 

Moreover, the obtained extended function (for which we keep the same nota- 
tion) can be extended by zero to the whole M d , outward the boundary e~ 1 ^d^l. 

Proposition 3.1. The spectrum of problem (|3.9p is real, discrete and consists 
of a countable set of points 

< a\ < a% < ■ ■ ■ < u e - < > +oo. 

The corresponding eigenf unctions can be normalized by 

W £ (v £ ,v £ ) = 5 l:i , (3.11) 

with W £ {u,v) defined by (pTTOj) . 

Proof. For any fixed e > 0, the bilinear form W £ (-, •) defines an equivalent scalar 
product in Hq (O e , e^^dQ). For brevity, we denote 

Hl w (Q £ ) = {wGH^(n £ ,e- 1 /' l dQ): \\v\\l w = W £ (w , w) < ^} . (3.12) 

Let G £ : L 2 (ft £ ) — > W {Q £ ) be the operator defined as follows: 

W £ (G £ f,w) = (f,w) L2{ ^ r wGH^e-^dQ). (3.13) 

Obviously, G £ is a positive, bounded (uniformly in e), self-adjoint operator. Since 
Hq W (Q £ ), for each fixed e, is compactly embedded into L 2 (Q), then G £ is com- 
pact as an operator from L 2 (Q) (Hq W (Q £ )) into itself. 

Thus, the spectrum o~(G £ ) is a countable set of points in R which does not 
have any accumulation points except for zero. Every nonzero eigenvalue has 
finite multiplicity. To complete the proof of the proposition, it is left to notice 
that in terms of the operator G £ the eigenvalue problem (|3.7j) takes the form 

G £ v £ = —v £ . 

□ 



l' £ 



We proceed with auxiliary technical results that will be useful in the sequel. 
Define the following norms in iif 1 (O e ): 



U>l| 2 



J a £ VvX7vdz-^ J Iv^dz + e 1 ^ J z)\v\ 2 da z ; 
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\\ v \\l^ = j a £ Vv -Vvdz + -^= J (x(e 1/4 z) - x(0)) \v\ 2 dz; (3.14) 

\\v\\ 2 s ^ Q = j a £ Vv -Vvdz + J (z T Qz)\v\ 2 dz. 

Lemma 3.3. The norms \\ ■ \\ e ,w> II ■ \\e,x an d II " \\eQ are equivalent. Morover, 

Ci \\v\\ 2 < \\v\\l w < C 2 \\v\\l „; 

(3.15) 

C 3 \\v\\l H < \\v\\l Q < C 4 \\v\\l x , 
with constants C\,Ci,C-& and C4 that do not depend on e. 
Proof. Indeed, by Lemma 14.21 and by the Poincare inequality, 



2 



\iJ <c- e V4 W ]i Vt ,]| ra ,~, Kd^Ml 



\ v \\e,W ~ \\V\\e,x\ ^ ^ II^H^^) || v v|| L 2 ( n;) 



and, thus, the first inequality in f|3. 15|) holds for sufficiently small e. 

The second inequality follows easily from the hypothesis (H4) and Lemma[ 

□ 

Remark 3.4. If v E ff 1 (O e ) decays exponentially, namely, 



\ v \\L 2 (R d \B R (0)) 



< Me" 7011 , 



for some constant M , then the norms defined in Lemma \3.3\ are asymptotically 
close. In particular, the following estimate holds: 

with the constant C = C(M, 70) independent of e. 

Lemma 3.4. Let /if be the first eigenvalue of the spectral problem (|3.7p . Then 
there exist two positive constants C\ and C2 such that 

Ci < iA < c 2 . 

Proof. The upper bound follows from (|3,8p and Lemma 13.11 The lower bound is 
the consequence of the boundedness of the operator G £ (see the proof of Propo- 
sition ED- □ 



3.2.1 Formal asymptotic expansion for the rescaled problem 

Following the classical asymptotic expansion method and bearing in mind Lemma[3 
we seek for a solution of problem (j3.7l) in the form of asymptotic series 

jjf = fjt + e 1 ^ + £ l/2 + . 

4 2 

V* = v ( z ) + £ l/4 Vl _ ( Zj Q + £ l/2 Vi _ ( Z; Q + £ 3/4 ^ ^ c) + 



(3.16) 
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where the functions Vk (z, £) are ^-periodic in £, k = 1, 2, ... . 

4 

Substituting ansatze (|3. 16[) into (|3.7p and collecting the terms of order e -5 / 4 
and e^ 1 in the equation, and of order the boundary condition, we 

see that the functions ui and vi do not depend on C. Then, collecting the terms 

4 2 

of order e _3//4 , we obtain that 

V3(>,C) = N k (()d k v(z) +w 3 (z), 

4 

where the vector function N(Q solves the problem 

' -div c (a(C)V c iV fc (C)) = div c a. fc (C), k = 1, d, C G V, 

aV^Nk ■ n = -dikUi, ( G S°, (3.17) 
[ AT fe (C) G ffi (y), 

The effective spectral problem comes out while collecting the terms of order e° 
and writing the compatibility condition for the resulting problem. It reads 

- div(a cff Vw) + (z T Qz) v = n v, v G L 2 (R d ), (3.18) 

where o cff is given by 

a|f = |7T / a ^(y)(^' + ^^i) ( 3 - 19 ) 
y 

The effective problem describes the eigenvalues and eigenfunctions of d-dimensional 
harmonic oscillator. In M 1 an explicit solution can be given in terms of Hermite 
polynomials. In the case under consideration we prove the following statement 
that characterizes the spectrum of problem (|3.18p . 

Lemma 3.5. The spectrum of the effective problem (|3.18p is real and discrete 

< [Ml < [12 < ■ ■ ■ < fij • • • — > +oo. 
The corresponding eigenfunctions Vj(z) can be normalized by 

(vi,Vj)g = J a cB Vvi ■ Vvj dz + j \z T Qz)viVj dz = 5ij. (3.20) 

R d R d 

We omit the proof of Lemma 13.51 which is classical. 

It is well known that the eigenfunctions of the harmonic oscillator operator 
have the form 

W'-, K ^' /2 f >-' /2 , (3.21) 

where Pk(z) is a polynomial of degree k. 

To summarize, the formal asymptotic expansion for v e takes the form 

z 



V 



where v is an eigenfunction of the limit spectral problem (I3.18p . is a periodic 
vector function solving (|3.17p . 

Notice that we can neglect the summands vi and ui since they do not depend 

4 2 

on the fast variable and thus, their iif 1 -norm is of order e 1 / 4 . 
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3.2.2 Justification 

Denote J(j) = min{i E Z + : fa = and let Kj be the multiplicity of the jth 
eigenvalue [ij of the harmonic oscillator operator (13.181) . 

The main goal of this section is to prove the following theorem. 

Theorem 3.1. Let hypotheses (H1)-(H4) be fulfilled. If(fip,v £ ) stands for pth 
eigenpair of problem \3. 7\ then the following statements hold true: 

1. For each j = 1,2, ... , there exist Ej > and a constant Cj such that the 
eigenvalue jx^ of problem (|3,7p satisfies the inequality 

Hj\<Cje l/4 , se(0,£j), 

where \ij is an eigenvalue of the harmonic oscillator operator (|3.18p . 

2. There exists a unitary Kj x Kj matrix (3 £ such that 
J00 +( %._1 



v p' 



k=ju) 



< Cj e x l\ p = J(j), • • • , J(j) + Kj - 1, (3.22) 



where 

ye = Vk{z) + £ 3/4 N ^ . Vvk{z y (3 33) 

Here the vector function N[Q solve problem (|3. 17)) : eigenf unctions Vk(z) 
of the limit problem are defined in (|3.18p .- the norm \\ ■ \\ £< q is defined just 
before Lemma [3~3\ 

Moreover, almost eigenf unctions {V?} satisfy the following orthogonality 
and normalization condition: 

J a £ Vn ■ VV£ dz + J (z T Qz) V k £ dz - 5 km \ < Ce 1 ' 4 . (3.24) 

Proof. The justification procedure will rely on Vishik's lemma about "almost 
eigenvalues and eigenfunctions" (see, for example, [9] and [5], p. 319, Lemma 1.5). 
For the reader's convenience, we formulate the mentioned result. 

Lemma 3.6. Given a self- adjoint operator IC £ : rl — > rl with a discrete spectrum, 
let v E K and v E rl be such that 

||u||% = 1, S = \\IC £ v — vv\\u < W\- 

Then there exists an eigenvalue fj,f of the operator K, 6 such that 

\nl - u\ < S. 

Moreover, for any 5\ E (5, \ v\) there exist {a £ } E M. such that 

where the sum is taken over the eigenvalues of the operator JC e on the segment 
[u — 6\, u + Si], and {uj} are the corresponding eigenfunctions. The coefficients 
a £ j are normalized so that ^ \ a j\ 2 = 1- 
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Let Hj be an eigenvalue of the effective problem (|3.18p of multiplicity Xj 
that is fj,j = fij+i = ■■■ = fj,j +Xj _i, and {v p (z)}, p = j,--- , j + Xj — 1, be the 
eigenfunctions corresponding to fij. Denote 

V £ (z) = v p (z) xe(z) + e 3/4 Xe(z) N{^) ■ Vv p (z), (3.25) 

where v p is the pth eigenfunction of the limit spectral problem (|3.18p . N is a 
solution of (|3.17p : Xs(z) is a cut-off which is equal to 1 if \z\ < £ ^ dist(0, dQ), 

—1/4 

equal to if \z\ > ^ — dist(0, dQ), and is such that 

< Xe(x) < 1, |V Xe | < Ce l l\ (3.26) 

We apply Lemma [3761 to the operator G £ : H^ W (Q £ ) — > Hq W (Q £ ) constructed 
in Proposition 13.11 (see (|3.13p ). The normalized functions V p = V p /\\V p \\ £j w and 
the numbers jij will play the roles of v 6 % and v 6 R in Lemma 13.61 Notice 
that Vj need not be equal to zero on the boundary e -1 / 4 <9f2; the cut-off function 
has been introduced in order to make approximate solution (|3.25p belong to the 
space Hl w (n £ ) (see (137121) ). 

Lemma 3.7. "Almost" eigenfunctions V £ are almost orthonormal. Namely, the 
following inequalities hold: 

(3.27) 

where W £ (u,v) and (-,-)q are defined by (|3.10p and (|3.20p . respectively. 
Proof We calculate first the gradient of the function V p . 

W p e = J £ lp {z) Xe(z) + e 3/4 J e 2p {z) + 4 p (z) Vxe(z), 

where 

Jf p (z) = Vv p (z) + V c (iV(C) • v p (z)) ; 

C=z/e J/ 

Jl p (z) = X e(z)V x (N(0-Vv p (z)); 
4 p {z) = v p {z)+e^N{^).Vv p (z). 

One can show that 



W £ (V p £ ,V £ )-J a £ (z)(xe(z)) 2 J £ lp -J[ 



lq dz 



Si: 



+^ Jv p (z)v q (z) (xe(z)) 2 dz 



<Ce l '\ 
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On the other hand, using Lemma 14.51 exponential decay of the eigenfunctions 
Vp(z) and the normalization condition (13.20p . we can prove that 



a £ (z)(xe{z)) 2 Jf p - J e lq dz 



x(0) 



V p (z)Vg(z) (Xe(z)) 2 dz 



+ £ V4 J q (e^) Vp { z ) Vq {z) (xe(z)) 2 da z - 5, 



pq 



<Ce l '\ 



Combining the last two estimates, we get 

\W £ (V p £ ,V £ )-5 pq \ <Ce l /\ 



which is the first estimate in (|3.27p . 

The second estimate in (|3.27p follows from the first one and Remark [3? 



□ 



Lemma 3.8. Let V £ = V £ /\\V £ \\ et w with V p defined by (|3.25|) , Then the follow- 
ing estimate holds: 

\\G e V £ p -{ H )- l V £ p\\ e , w <Cpe l l\ p = i,— ,% + Kj-l. (3.28) 

Proof. Simple transformations result in the following relations: 



\\g £ v £ - ^i)- 1 v;\\ £ , w = \\vf\\-ly sup w £ {g £ v £ - (^r 1 v £ , 



w 



W e.W = l 



By (M), 



1 



\\g £ v £ p - {^r l v;\\ e , w = f ||y p £ ||-^ sup {^ p (v;,w) L2m 



Vp \\w\\ SiW =l 



a £ W p -Vwdz + 



^ jv £ wdz- e 1 / 4 J q(eV*z) V £ w da z } 



-\\Vp e \\^ w sup {/? + If + £ 3 / 4 If}. 
Vp 



\\w\L.w=l 



Here 



I t = Vp j Xe{z) Vp(z) w(z) dz - J (x(e 1/4 z) - x(0)) Xe(z) v p {z) w(z) dz 



(C) (Vvp(z) + V c (iV(C) • Vvp(z))) ■ Vw Xe(z) | c=z/e3/4 dz; 
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If = -i= y x(e 1/4 z) f p (z) x e («) cfa - e 1/4 ^ g(e 1/4 2:) u p (*) w(z) da; 



i I = Vp I Xs{z)N{-^ 1 ) ■Vv p (z)w(z)dz 



z) ■ VwUp(z)| c=z/e3/4 dz 



- j a(()V z (xe(z) N(C) ■ Vv p (z)) ■ Vw 



C= 2 / e 3/4 dz 



+e 1 / 4 x(0) j X£ (z)N(^ JI ).Vv p (z)w(z)dz 



-e 1 / 4 J q(e 1 / i z) X£ (z)N(-^) .Vv p (z)wda. 

Integrating by parts in the last integral in If, taking into account (H4), (|3.17p 
and (|3.26p . we obtain 



+ J div,(a(C)(/ + V ( N(())Vv p (z)) \ Q=z/em w(z) X e(z) dz + O^ 1 / 4 ), e -> 0. 

Here we have also used Lemma [3731 and the fact that ||u;|| Sj vf = 1. 

Bearing in mind the definition of the effective diffusion (|3.19p and (I3.18p . by 
virtue of Lemma 14.51 one has 



\n\<CeV A \\w\\ H1{Rd y (3.29) 

By Lemma 14.21 



\I £ 2 \ < Cae 1 / 4 \\v p \\ H i m \\w\\ HHRd) . (3.30) 

Using the boundedness of and the regularity properties of N,v p ,Xe, one can 
show that 

|If| < C 3 \\Vv p \\ L 2 (Rd) ||iw||iri(R*). (3.31) 
Using Lemma 13.71 we see that for small enough e, 

\K% W > \ (3.32) 

Finally, combining (|3.29p - (|3.32p we obtain the desired estimate (|3.28p . Lemma [3.8l 
is proved. □ 
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By Lemma 13.61 in view of the estimate obtained in Lemma 13.81 for any 
eigenvalue fij of the effective problem (|3.18j) there exists an eigenvalue of the 
original problem such that 

K (i) -^|<C^ 1/4 , (3.33) 

where q(j) might depend on e. 

Moreover, letting 5\ in the statement of Lemma 13.61 be equal to Qj e 1 / 4 (the 
constant Qj will be chosen below), we conclude that there exists a Kj(e) x kj 
constant matrix a £ such that 

Ji+Kj(e)-\ C£l/A 



n- 3 E 



Y kp u k 



< 2—— < CjiQj)- 1 P = j,--- J+nj-l, (3.34) 

e,W d l 



here fi £ k , k = Jj(e), . . . , Jj{e) + Kj(e) — 1, are all the eigenvalues of operator 
(Cr £ ) -1 which satisfy the estimate 

|/4-Mil<e ie V4 (3.35) 

Since the eigenvalues fij do not depend on e, one can choose constants Sj > so 
that the intervals (fij — QjE 1 ^ , fij + Qj£ 1 / 4 ) and (/ij — ©je 1 / 4 , fn + ©jE 1 / 4 ) do 
not intersect if fij / /x, and e < mm(ej, Si). Then the sets of eigenvalues {fi%} 
related to different fij in (|3.35p do not intersect for sufficiently small e. 
In the following statement we prove that Kj(e) > Kj. 

Lemma 3.9. The columns of the matrix a £ , that is the vectors {oi £ p} p=J ^ J 
of length Kj(e) are linearly independent. As a consequence, Kj(e) > Kj. 

Proof. A simple transformation gives 

Jj+Kjie)-! 

w £ (v £ , v £ q ) = w £ (v £ - E <V! -vm 



Jj+Kj(e)-1 Jj+Kj{e)-l Jj+Kj(e)~l 



+w £ ( E a% p v £ f,vi- E + E «I P <v 

Taking estimates (|3.27p and (|3.34p into account, we obtain 

J, +Kj(e)-1 

E a kp a kq-$p,q ^C&J 1 , p,q = J(j),--- ,J(j) + Kj -1, 
k=Jj 

and, in other words, 

\(a £ p ) T a £ - 5 p , q \ < CQJ 1 , p,q = J(j),--- ,J(J) + Kj - 1, (3.36) 



where a £ p denotes a pth column in the matrix a £ . The last inequality means 
that the vectors {a £ p } p ^j^ J 1 are asymptotically orthonormal, as Qj grows 
to infinity. This property implies the linear independence of the vectors {a £ p } 
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for sufficiently large Qj. Indeed, assume that {a £ p }3oo J 1 

are not linearly 

independent. Then there exist constants Cj^, • • • , Cj^ +K3 _i such that 

j(j) +Kj -i 

Without loss of generality we assume that Cj^ = 1 > max*, \ci-\- Then 



a £ JU \ + ^2 Ck a% = 0. 



k>j(j) 

Multiplying the last equality by a £ w.n and using (|3.36p we obtain the inequality 

that contradicts (|3.36p if Qj 1 is sufficiently small. Thus, the vectors {a^j^j^ 1 
of length Kj{e) are linearly independent. Obviously, it is possible only in the 
case Kj(e) > Kj. □ 

Lemma 3.10. For any q, < m < fi e < M q . 

Proof. The estimate from below is the immediate consequence of the bounded- 
ness of the operator G e constructed in Proposition 13.11 

To obtain an upper bound for [i £ q , we recall estimate (|3.33p . For any j, there 
exists an eigenvalue of problem (13.71) converging to the jth eigenvalue of the 
effective problem. Namely, the estimate holds 

K*(i) _ ^J(J')I - C i £l/4 ' 
where J(j) = min{z 6 Z + : [ii = fij}. Obviously, q £ {j) > J(j)- Thus, 

that implies the desired bound. □ 

Our next goal is to prove that any accumulation point of the sequence /j, q , as 
e — > 0, is an eigenvalue of (I3,18p , 

Lemma 3.11. If, up to a subsequence, [i e q — > fi* , as e — > 0, then fi* is an 
eigenvalue of the effective spectral problem (|3,18p . 

Proof. Since fj, q is bounded, then 

lbg||e,VP < c q 

with || • \ \e,w defined in (|3.14p . In view of Lemmata 13.31 and 14.41 the eigenfunction 
v q (extended to the whole M. d ) converges weakly in H 1 (M. d ) and strongly in L 2 (M. d ) 
to some function v* . To prove that (fi* , v*) is an eigenpair of the effective problem 
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(|3.18|) . we pass to the limit in the integral identity (|3.9p . Using standard two-scale 
convergence arguments we obtain 

J a cB Vv* -Vwdz + J (z T Qz) v*wdz = // J v* w dz, w £ H 1 ^). 

R d R d R d 

The last equality is the weak formulation of ()3.18j) . Since \i e — > /j,*, as e — > 0, 

then considering (13.911 and (13.111) we conclude that lim lluf ||? ,p; , = a* . Using 
* ' y e->o" q "L 2 (n e ) 



the strong convergence of v 1 ^ in L 2 (M. d ), we see that ||f*||^2( K d) > M*- By Lemma 
we have fx* > 0. Therefore, v* ^ 0. This completes the proof. □ 



Lemma 3.12. Let fij be the jth eigenvalue of problem (I3.18P of multiplicity Kj, 
that is jjLj = fij+x = ■ ■ ■ = fij +Kj _i. Then there exist exactly Kj eiegnvalues of the 
original problem (|2.ip converging to it. 



Proof First, we prove that there are not more than Kj eigenvalues of problem 
(j3.7|) converging to fj,j. Assume that there exist Kj + 1 eigenvalues fJ-j^j) such 
that 

VJ^ti) ~* H> fc = 1, • • • , Kj + 1. 

By Lemma |3.11| the corresponding eigenfunctions Vje^, extended to the whole 
converge weakly in i? 1 (R d ) and strongly in L 2 (K d ) to the eigenfunctions vt 
of the effective problem (|3.18p . k = 1, • • • , + 1. Passing to the limit in the 
normalization condition (|3.1ip yields 

«> u fc)L2(R<<) = r^r-lkSik, i,k = 1, - ■ ■ ,Kj + l. 

I * \d 

Therefore, eigenfunctions {vt}i.Li corresponding to fij are linearly independent. 
Recalling that the multiplicity of fij is Kj, we arrive at contradiction. Thus, there 
are not more than Kj eigenvalues of problem (|3.7p converging to \ij. 

On the other hand, by Lemma [3791 there exist at least Kj eigenvalues of (|3.7p 
converging to fij of multiplicity Kj. Lemma 13.121 is proved. □ 

Combining Lemmata I3.101f37l2l completes the proof of the first statement of 
Theorem 13.11 

We turn to the proof of the second statement in Theorem 13. 11 

First of all, let us notice that the orthogonality and normalization condition 

(|3.24p follows directly from Lemma 13.71 and the exponential decay of Vk(z) as 

eigenfunctions of the harmonic oscillator. 

In order to prove estimate (I3.22p . we recall the estimate obtained in Lemma l3.8l 

and apply the estimate in Lemma 13.61 with 5\ = Cj, Cj being a sufficiently small 

constant. This estimate reads 



C e i/4 

<2-— ^e 1 / 4 P = j;--- J + Kj-1, 
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where S(j,e) is the set of eigenvalues fii satisfying the estimate 

1/4 — Hj\ < Cj) 

the constant matrix a £ is such that 

\(a £ p ) T a £ q - 5 Ptq \ < Cj e 1 / 4 , p, q = J(j), • • • , J(j) + Kj - 1, (3.37) 

Prom the first statement of Theorem 13. II we deduce that the set S(j,e) coincides 
with the set of eigenvalues \ f° r sufficiently small e. Therefore, 



E 

k=j(j) 



< Cj£ 1/4 p = j;--- ,j + Kj - 1, 



(3.38) 



e,W 



with a constant Kj x Kj matrix a 6 which satisfies inequality ()3.37p . 
It remains to use the following simple statement. 

Lemma 3.13. For any n x n matrix A satisfying an equality 

p T A-I|| £(R „ iR n } = 7 G (0,1), 
there exists a unitary matrix B such that 

\\AB - V\c(m. n ,Mp) < 7; 

here I is a unit matrix, and 

\\D\\c(r^,r^) = sup \\D£\\. 

llfll=l 

We omit the proof of this lemma which can be found in [S]. According to 
()3.37p and Lemma |3.13| there exists a unitary Kj x Kj matrix (3 E such that 

WofF-nc^^^Cje 1 /*. (3.39) 

Taking into account Lemma 13.71 estimates f|3.38|) . (|3.39p . one can show that 

j(j)+ Kj -i 

\- Yl ^P V k <C 3 e l '\ p = J(J),-.. ,j(j) + Kj -i. 

k=J(j) 

Due to the exponential decay of the eigenfunctions Vk{z) defined in (|3,18p . one 
can replace VE defined by (|3.25|) with (|3.23p . Then, by Lemma 13.31 similar 
estimate holds for || • \\ e q norm. Theorem 13.11 is proved. □ 



Bearing in mind the result obtained in Theorem 13.11 we formulate the main 
result of the present paper characterizing the asymptotic behaviour of eigenpairs 
(Af.uJ) of problem (l2~Tp . 
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Theorem 3.2. Let conditions (H1)-(H4) be fulfilled. If {X £ -,u £ -) stands for the 
jth eigenpair of problem (12.ip . then for any j, the following representation takes 
place: 



d-l 



q(0) + ^= £ , ^(x)=^(- T7i ) 



where the eigenpairs (fj, £ - , Vj (z)) of problem (|3.T|) are such that 

1. For each j = 1,2,..., there exist Ej > and a constant Cj such that 

— Hj\ < Cj e 1 / 4 , e G (0, ej), 

where fij is an eigenvalue of the harmonic oscillator operator ()3.18j) . 

2. Let fij be an eigenvalue of (|3.18p of multiplicity Kj, that is \ij = ■ 
fj,j +Kj _i. Then, there exists a unitary Kj x Kj matrix f3 £ such that 



fc=j(j) 



^•e 1 / 4 , p = J(j),-.. ,J(j) + «,-!, 



where 



V k £ = v k (z) + e 3 / 4 N{-^).Vv k (z) 



Here the vector function N(Q solve problem (|3.17p ; eigenf unctions Vk(z) 
of the limit problem are defined in (|3.18p ; f/ie norm \\ ■ \\ £ q is defined in 

BP- 



4 Auxiliary results 

Lemma 4.1. For any u> e (x) E -^(f^, <9f2) t/ie following estimate holds 



1 IS 



d-l 



e y 



Iu7 e | 2 dcr 



<C\\w £ \\ L 2 {ne) \\Vw £ \\ L 2 {yis) 



wii/i a constant C independent of e. 

Proof. Introduce a y-periodic vector function x{v) as a solution of the following 
problem on the periodicity cell Y: 



-divyX 



IS 



d-i 



\y\d 

( X ,n) = -1, yeS 

Notice that \ is a smooth function. Then 

i . (X, |S°U-i 
-^iv,x(-)=^ 

„s|2 



y ey 



Multiplying the last equality by \w £ \ and integrating by parts over Q £ yields 



1 IS 



w £ \ 2 dx 



\w £ \ 2 da 



. Ivl , • ! , • , , (x{-),V\wr)dx 

£ I * \d J J J 

iic ^£ S2e 



that easily implies the desired estimate. Lemma is proved. 



□ 
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Lemma 4.2. Let Q £ = e~ a Q £ , T, £ = e~ a T, £ . Then, for ip(z) G H%(Q e , e" Q 3ft) 
and tp G C 1 (]R rf ) J t/ie following estimate holds 



1 IS 



d-l 



V{s a z)\il){z)\ 2 dx- / ^(e Q z)|V(z)| 2 da 



e i-« |y| d 
^cilV'll^^llWll^^)- 

wit/i some constant C independent of e. 

Lemma 14.21 is proved in the same way as Lemma 14.11 

Lemma 4.3. Suppose two nonnegative functions /i,/2 G C S (B), defined on a 
bounded domain B, are such that x = is the global minimum point for both of 
them, and /i(0) = /a(0) = 0. Moreover, assume that 

H(f k )(0)>al, a>0, 

where H{fk) is the Hessian matrix of fk, k = 1, 2. 
T/ien i/iere exists a constant C such that 

c h < h < c- 1 fi. 

Proof. Assume that there exists a sequence Xj G B such that 

fl(%i) 

iiyi _> o, j oo. 

Since /2 is bounded, then f\{xj) — > 0, as j — )• oo. And, consequently, Xj — > 
0, as j — » oo. Recalling that H(f i)(0) is bounded from below, we arrive at 
contradiction. Lemma is proved. □ 

Lemma 4.4. Compactness result 

Denote 

H^(M. d ) = {u; G fl' 1 (M d ) : \\vf Q = J \Vv\ 2 dz + J (z T Qz)\v\ 2 dz < 00}. 

Then i^g(]R d ) is compactly imbedded into L 2 (W 1 ). in oi/zer words, any {v n } C 
^(R^) suc/i that \\v n \\Q < C, converges strongly along a subsequence in L 2 (R d ). 

Proof. Obviously, v n , up to a subsequence, converges weakly in L 2 (R d ) to some 
function v* , n — )• 00. Let us prove that 11^11^2^) — > \\v*\\i?md), as n — > 00. 
Since 

\z T Qz)\v n \ 2 dz < C, 



one can show that for any 5 > 0, there exists a ball B R ^(0) such that 



l d \-Bi* W (0) 



| 2 cte < 5. 
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Without loss of generality we assume that ||un||z, 2 (R d ) = 1- Then 

\\ V n\\h(B R(s) (0)) = 1 " W V n\\ 2 L2(m<l)\B R(s) (0)) > 1 ~ & • (4-1) 

Since IKIIin-CB^o)) < C> then IK ~ v * \\l 2 (B r(S) (0)) -> 0, as n -)• oo. Passing 
to the limit in (|4.ip . we have 



|L2 (R d) > |k*|| L 2 (Bjiw( 0)) 



On the other hand, 



\ v \\L 2 (R d ) 



< liminf ||iy, 

n— >oo 



> 1-5 2 . 



1. 



Combining the last two inequalities yields \\v* Wipngd^ = 1. Lemma is proved. □ 
Lemma 4.5. Mean-value theorem 

Let <3? G £ 2 (y) be sitc/i that f Y $dy = 0, and V G C 1 ^) sate/?/ i/ie estimate 

\D k V{z)\ < Ce- 7o|z|2 , 7 > 0, fe = 0,1- (4.2) 

Denote by xi x ) a cut-off which is equal to 1 if \x\ < gdist(0, dti), equal to if 
\x\ > f; dist(0, dfl), and is such that 

< X (x) < 1, |Vx| < C. (4.3) 

Then the following estimate holds: 

I ^(-J V (^)x(x)W(^)dx ^Ce^e^Wn^WWy^ 

J c c, c 



for any W G i? 1 ! 

Proof. Since f Y & dy = 0, then there exists a periodic vector function <p(y) such 
that 

-div„^(y) = y G y 

(^,n)=0, 

and |M|_L2(y) < C || $|| £,2 (yy Changing the variables we have 

-£divy>(~) =*(j). 

Multiplying the last equation by V(p?) XeO^) W(^r) , integrating by parts over 
O e and using (fOD we get 

*(f)^)x(x)W(^)cfa 



V(f)-V[^)x(x)w(^) 



r/.r 



< Ce l-a £ da /" |^_i_)| e -70M 2 [|^| + dz 



n e 



< Ce L - a e 



l—a Aa 



\\W\\h^)( H^)\ 2 e- 2 ^ 2 dz) 



1/2 



0= 



< 



Cs^e^WWy^yh^yJ / e~ 2 ^ dz 



1/2 
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The integral in the parenthesizes converges. Lemma 14.51 is proved. 



□ 
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